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based X-ray telescopes. We consider axinos and gravitinos with masses between 1 keV and 100 
keV in supersymmetric models with small R-parity violation. We show that axino dark matter 
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end of the allowed range of the axino decay constant f a , while the gravitino dark matter remains 
unconstrained. 
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I. INTRODUCTION 



Though many popular supersymmetric (SUSY) models tout the lightest neutralino as 
the dark matter candidate, the axino or the gravitino, too, may qualify for a similar role. 
The stability of these dark matter candidates is generally ensured by the conservation of im- 
parity, defined by R = (— i)( 3B + L + 25 ) ; where B, L and S are baryon number, lepton number 
and spin respectively i?-parity keeps the proton from decaying. As a by-product, it keeps 
the lightest supersymmetric particle (LSP) stable. 

This, however, need not always be the case. It is sufficient for the dark matter particle 
to be stable on cosmological scales rather than be fundamentally stable. Thus a slowly 
decaying light axino or gravitino in an R-parity violating theory is not incompatible with a 
universe containing dark matter. Such a situation has a distinct phenomenology of its own. 

Since the simultaneous violation of B and L would lead to unacceptably fast proton decay, 
one usually assumes only one of them to be violated by odd units. L- violation, in particular, 
leads to useful mechanisms for the generation of neutrino Majorana masses at both the tree- 
and one-loop levels. 

Prima facie, L-violation by odd units renders the LSP unstable. This rules out the 
candidature of the lightest neutralino as a dark matter constituent as it tends to decay fast 
compared to the cosmological scales. On the other hand, if the LSP is a light gravitino 
or an axino, it is still possible to have a viable dark matter candidate. The requirement 
is that the decay of the gravitino or the axino through R-parity violating couplings should 
be sufficiently suppressed such that their lifetime is larger than the age of the universe. In 
this note, we explore this possibility and study the consequent astrophysical observations, 
in terms of X-rays coming from galaxies and galaxy-clusters. 

The axion field in a supersymmetric model is the pseudoscalar part of an electrically 
neutral chiral supermultiplet $ a . $ a can be expanded as $ a = (s + ia)/y/2 + 6a + 99F, where 
the scalar field is denoted by s, the axion, or the pseudoscalar field, by a and their fermionic 
partner, the axino, is denoted by a. As long as SUSY is unbroken, the supermultiplet 
remains light, since it is protected by the Peccei-Quinn synmetry {7(1) pq l|-|3j. The masses 
split, lifting the axino from the almost massless axion, if SUSY is broken. This is achieved 
either at the tree level or via loop diagrams. The precise value of the axino mass depends 
on the model. While one would typically expect the axino mass to be around Msusy, the 
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soft SUSY breaking scale, it has been realised that it can be much smaller than 

-^susy depending on the mechanism of mass generation. In the present work, we set the 
axino mass as a free parameter, which could span many classes of models. 

It is interesting to note that the axino is a good dark matter candidate even if R-parity is 
not exactly conserved. The lifetime of the axino can be very long because of the suppression 
from a large Peccei-Quinn (PQ) symmetry breaking scale and a small amount of R-parity 
violation. Axino cold dark matter (CDM) with R-parity breaking has been considered in 
[5I and subsequently in the context of the INTEGRAL anomaly in 6]. For non-thermal 
production, the axinos with a mass in the range 2-100 keV can be the warm dark matter 
(WDM) 0,3- 

Another candidate for warm dark matter in supersymmetric theories is the gravitino. Al- 
though the mass of the gravitino in the simple types of supergravity (mSUGRA) models is 
in the electroweak/TeV scale, situations with a very light gravitino, too, are often envisioned 
in supersymmetry. For example, light gravitinos are natural in models with gauge mediated 
supersymmetry breaking. Such a light gravitino can be the dark matter candidate under 
certain circumstances 9]. These light gravitinos can lead to novel signatures jl], lo| at col- 
iders and accelerators and further have strong astrophysical and cosmological implications 

The light LSP's we are considering i.e. the gravitino /axino could form either cold dark 
matter or warm dark matter depending on their mass ranges. Typically a ~GeV dark matter 
would be a CDM whereas a keV dark matter would be a WDM. The exact limits however 
depend on the details of the candidate, its couplings etc. An axino with masses between 2 
keV to 100 keV has been shown to be consistent as a warm dark matter candidate satisfying 
the relic density constraints |3] • In the present work we will consider this mass range for the 
axino, however, with R-parity violation. 

A particularly simple and minimal model of R-parity violation is the case where there 
are only bilinear AL = 1 terms in the superpotential, 

Wf = eiLiH u , (1) 

where L, and H u are the leptonic and up-type Higgs super-fields and the index % = {1, 2, 3} 
runs over the generations. The model which has been studied over the years from many 



angles 11] is specified by six additional parameters at the weak scale: the three sneutrino 



vacuum expectation values (VEV) and three bilinear R-parity violating couplings. The 



presence of these bi 
the width given as 



inear terms leads the axino to decay into a photon and a neutrino with 



7] 



1 TH ^ 

r(S -> 1Vl) = » ^ jf «™ & ( 2 ) 



where = (v^/v with the Higgs VEV given by v = 174 GeV, 



N- N* 



AO 



is the photino-zino mixing parameter containing the neutralino mixing matrix N and the 
mass eigenvalues m^, a em is the fine structure constant, Mz is the mass of the Z-boson, 
rria is the axino mass, f a is the Peccei-Quinn symmetry breaking scale (up to a domain wall 
number), also identified with the axion decay constant, and C ail is axion-photon-photon 



coupling, a model dependent quantity (See Kim |12| for a table of values). If we restrict 
ourselves to hadronic axion models, the axino decay into a photon and a neutrino arises 
from the anomaly coupling with the photon vector multiplet and the neutralino-neutrino 
mixing generated by the sneutrino VEV (z/j) Q. 

For a suitably small mixing parameter, it is remarkable that the lifetime of the axino 
can exceed that of the age of the universe, for a Peccei-Quinn symmetry breaking scale of 
(^(lO 11 ) GeV and an axino mass of 5 keV. This can be seen from the life time expression for 
the axino given as: 



&\ U jz\\ 2 / rn h \ -3 / f a 



In a similar manner, gravitino can also live long enough to be considered as dark matter 
though it is governed by a different set of couplings. In a recent analysis by Buchmueller 
et. al [isl it has been shown that gravitinos of mass (9(5 GeV) or just above are compatible 
with thermal leptogenesis, cold dark matter and Big Bang Nucleosynthesis. However, purely 



from a dark matter relic density point of view, it has been shown by Rubakov et. al [14 1 
that a lighter gravitino of the order of a keV could as well be consistent with relic density 
observations as a warm dark matter candidate. In the presence of bilinear R-parity violation, 



gravitino decays into a neutrino and a photon 



13 



expression [16[ : 



151 ] . The decay rate is given by the 



^"►•H^l^l*^, (5) 



where m 3 / 2 is the gravitino mass and Mp is the Planck scale, and U^ u is given by 13] 



10 8 (lO- 7 ) (200Gev) ~ 10 5 (lO" 4 ) (200Gev) 



(6) 

Here, g z is the gauge coupling, c^ v and c za are the mixing elements in the neutralino-neutrino 
mixing matrix, v v is the sneutrino VEV, m is a common mass scale for the superpartners 
and €3 is the bilinear R-parity violating parameter defined in Eq.(l). These parameters 
appearing in the second term of the RHS can be found in 
lifetime of a 5 keV gravitino we find , 

2-body„, vln W ^3 rV ™ V ( 



13J. Using this formula for the 



T 3/2 



From the above we see that the lifetime of the gravitino is several orders larger compared 
to the age of the universe (C(10 18 s)). 

Since, according to the above expressions, their lifetimes are large but finite, one would 
expect a few of the axinos/gravitinos to decay through the two body processes (|2j) and (jSJ). 
The resultant photons are in the X-ray regime and thus could be observed by the astro- 
physical X-ray telescopes. These X-rays would constitute an excess over the astrophysical 
background. In the recent years, there has been a report of such an observation by the 
Chandra satellite based observatory with a line emission at 2.5 ± 0.11 keV (90% CL) and 
flux of (3.53 ± 2.77) x 10^ 6 photons cm" 2 s~ l at (90% CL) 0. 

Recently, Loewenstein and Kusenko 18] have tried to explain the Chandra excess [Tj 
through a radiative decay of a sterile neutrino of mass of 5 keV. It should be noted that 
astrophysical X-rays have been observed earlier too from the Suzaku satellite and have led 
to stringent constraints on sterile neutrino mixing and mass parameters In the present 
work, we would like to obtain similar constraints on the axino masses and the Peccei-Quinn 
symmetry breaking scale f a for a fixed value of the neutralino-neutrino mixing parameter 
satisfying neutrino data, using the Chandra results. On the other hand, similar analysis 
on gravitino dark matter does not provide any constraint because of the larger suppression 
through l/M Pl in the decay width. 
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II. DARK MATTER HALOS AND FLUX OF X-RAY PHOTONS 



The X-ray flux from the dark matter halos can be calculated for a given galaxy knowing 
its luminosity distance and the dark matter distribution. In the following we outline the 
calculation of the flux in general and apply it for our own galaxy. 

An object such as a field galaxy or dwarf galaxy or a cluster of galaxies possessing a dark 
matter halo of mass Mom will be composed of N = Mdm/^x dark matter particles of rest 
mass mx- If Tx is the dark matter particle decay rate into photons of energy E 11 then the 



total associated X-ray luminosity is 20] 



£ = ^LM DM T X . (8) 
m x 

Here we have assumed that the halo is relatively nearby and that redshift effects on the 
luminosity are negligible. 

We will also consider the simple case where E 1 = mx/2 (when the neutrino mass is 
negligible compared to that of the dark matter candidate), which is relevant for the two 
body decays which are the focus of the present work. The flux from an object seen through 
the telescope is simply 

£ 



47r£>!, 



(9) 



where Dl is the luminosity distance to the object. This can be taken to be the coordinate 
distance r from the earth to the dark matter site (neglecting redshift effect). 
With E 1 = mx/2 we get 

£=^X (10) 

Hence, the flux F can be written as 

= IT / \pDM(r)r 2 drdn, (11) 

Jl.o.s r 

where Pdm(^) is the density of dark matter particle as a function of distance from the 
galactic centre. In order to calculate the flux which is coming from the direction of the 
galactic anti-center, we have 



S7T 



/0dm (r) dr. (12) 

re 



Field of View (fov) 




FIG. 1: Solid angle (fi) calculation of the dark matter halo. 



The flux which is coming from the direction of the galactic center is given by 



S7T 



Pnu{r)dr. (13) 



Here r (~ 8 kpc) is the distance of an observer on the earth from the galactic centre (of 
Milky Way). Note that this derivation applies to integration along the Milky Way only. The 
reason we are considering Milky Way is because of the fact that it is a very extended source 
of the dark matter halo and SPI telescope (which has a wide field of view) is well suited for 
such an extended object. Our analysis on SPI telescope will be presented in Sec. IV and we 
shall show that this will give a very strong constraint on the axino dark matter model. 

In order to determine the dark matter halo contribution to the X-ray flux, one needs to 
know the distribution of the DM in the halo of the Milky Way. In Eq. (fTT|) . Q is the small 
(^C 1) solid angle corresponding to the field of view of the telescope which can be calculated 
from figure [TJ In figured], the projected radius of the field of view is denoted by Rf ov , and 
the angular field of view by 9. Approximating Q as the solid angle of a cone with an apex 
angle 9, we have 

n-£, d4) 



where 

9 = ^1^. (15) 



2R f or 



r 

-3 



In particular for ACIS aboard Chandra, 9 fa 5 x 10 rad, which determines the solid angle. 
Later on we shall show that a much larger solid angle is obtained for SPI. 

The final input needed in the calculation of the flux is Pdm ( r ) j the density profile of the 
dark matter. While the actual profile of the dark matter in the galaxy is unknown, there are 
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many profiles put forward in the literature. Many of them also have analytical expressions. 



For example, the isothermal profile 21] is given by 



PdmW = Vh 2 ] 2 , (16) 



where Vh corresponds to the contribution of the DM halo into the galactic rotation curve 
in its flat part, Vh ~ 170 km/s. Here the core radius r c is almost 4 kpc. Now, the profile 
integral in Eq.( fl2l) for the galactic anti-centre is 



\ pDM{r)dr= — 

Jr Q r c 



^-tan- ^ 
2 V r r 



^tan" 1 ^), (17) 



where C = v^/AiiGn- Similarly for the galactic centre case from Eq.([T3i): 



r° / x , C \tx /r G \ 



(18) 



Using Eqs.f ll7p and (118p one can calculate the flux from the galactic anti-center and galactic 
center, respectively. 

Once again let us mention that the calculation of the dark matter line flux in the direction 
of the Galactic center or anti-center has been performed just for simplicity. In general one can 
calculate the flux coming from any given direction on the sky. We shall present our analysis 
for an angular distance 13° between the given direction on the sky and the direction towards 
the Galactic center when we present our results for SPI. 

Another point is that in Eq.(12) and (13) the upper (lower) limit of integration is actually 
not infinte but depends on the assumed size of the Milky Way dark matter halo. However, 
since the integral is converging, we shall obtain the same result as long as the limit is 200 
kpc or above. This is the reason we have assumed the upper (lower) limit to be infinte just 
for simplicity. 



The data from Virgo cluster by ACIS aboard Chandra has been analyzed in Ref. [20 1. 
The instrumental background has been estimated to be 2 x 10~ 2 cts s" 1 in a 200 eV energy 
bin. It has been shown that in order to detect a dark matter decay line signal of energy 
Ej = mj(/2 at 4<r level, one must observe a flux above the value, Fdet = 10 -13 erg cm _2 sec _1 
with an integration time of 36 ksec for observation. In our analysis of the Milky Way dark 
matter halo, we shall assume that the same flux limit is applicable in the entire x-ray energy 
range of 1 keV to 8 keV relevant for Chandra telescope. This determines the region of the 
parameter space for any theoretical scenario, which can be probed by Chandra. 

9 




FIG. 2: Flux of photons from the Milky Way galactic centre for the Max. disk model The 
horizontal line is the lower limit of detectability of the flux applicable for the Chandra X-ray 
satellite observations as explained in the text. 



III. NUMERICAL ANALYSIS FOR CHANDRA TELESCOPE 



We have taken the value of the neutrino-neutralino mixing parameter (dU^z = U^ u ) to 
be 9 x 10~ 7 , which is the allowed value from neutrino data. In figure El we show the flux 
from the Milky Way galactic centre as a function of the axino mass m 5 , following the Max. 



disk model A 2 



22] of NFW profile 23|, |24| for various values of f a . As can be seen from 



the figure, the lower is the values of f a , the larger is the flux. From the figure we see that 
for the astrophysically unconstrained range of the Peccei-Quinn symmetry breaking scale, 



namely, f a > 10 9 GeV 



251 ]. the axino mass is unrestricted in the range shown in the figure. 



However, for lower values of f a some constraint can be obtained. 

The question that arises is: for which profiles there exist constraints from Chandra within 
the astrophysically unconstrained values of f a . In figure [3] we present the fluxes for various 
orofiles including the isothermal profile discussed in the previous section. The isothermal 



2l| , and the NFW profiles 22|, |26] do not provide constraint if one takes the lower limit on 
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FIG. 3: Flux of photons from the Milky Way galactic centre for various 



profiles. Top left is for 



221 ] . Similarly, bottom left 



isothermal profile, top right is for the NFW favored models (A\ or B\ 
is for the NFW Max. disk models B2 of the MW DM halo whereas bottom right is for the NFW 
profile best fit parameters given by Battaglia et al 26(. As before the horizontal line is the lower 



limit of the flux applicable for the Chandra X-ray satellite observations. 



fa ^ 10 9 GeV strictly. However, there are ways of avoiding these astrophysical constraints in 



some situations such as (a) existence of paraphotons |27l.l28j] and (b) existence of temperature 
dependent couplings 29|. In either case, f a lower limit can be relaxed, sometimes even 
to values of (9(10 6 ) GeV. In such situations the constraints from Chandra could play an 
important role. 

Before concluding let us briefly touch the constraints on our model coming from the 



observation of the Virgo clusters 



20 



30 



3l|. Using Eqs.g]), (|5j| and (fTU|), the flux from a 
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2 4 6 8 10 12 14 16 

m s (KeV) 

FIG. 4: Flux of photons from the centre of the Virgo cluster. As before the horizontal line is the 
lower limit of the flux applicable from the Chandra X-ray satellite observations. 



cluster can be written as, 

* - , 19 x ^ ergs cm -v (^) 2 ( (^) 2 ( 2 (19 , 

The values for Dl and Mqm f° r t ne Virgo cluster are given in [2(|. We have used Mgfa = 
1O 13 M and D^ = 20.7 Mpc in Eq. ffT9]) and demanded that the flux should be less than the 
detectable limit 10~ 13 ergs cm -2 s _1 . The results are plotted in figure HJ It is evident that 
the constraints obtained, are at least as stringent as those for the Milky Way, for the more 
optimistic density profiles. Since the luminosity distance -Dl is much larger for the Virgo 
cluster, one would expect that the corresponding flux will be smaller than that obtained for 
Milky Way. For the Virgo cluster, however, the total dark matter mass in the field of view 
(M^j) is much larger compared to that of the Milky Way in a particular direction (say, the 
direction of the galactic center or anti-center). This large mass for Virgo cluster compensates 
for the large luminosity distance and hence the resulting flux becomes comparable to the 
flux from Milky Way for more optimistic density profiles such as the NFW favored model 



A x or B 1 



22]. 
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IV. NUMERICAL ANALYSIS FOR SPI 



In this section we shall look at the bounds on the axino parameters originating from 
an analysis based on the data from the high-resolution spectrometer SPI on board the 
INTEGRAL sattelite 32|. SPI is sensitive to X-ray photons in the energy range 20 keV to 
7 MeV. Obviously, the data from SPI is applicable to axinos in the mass range of 40 keV 
to 14 MeV. However, we shall restrict ourselves to the axino mass range of 40 keV to 400 
keV, which includes the window on the axino mass where the axino can act as a warm dark 
matter candidate. The negative search of a dark matter decay line using SPI put restrictions 



321 ] . Using this flux restrictions we derive bound on the 



on the dark matter decay line flux 
axino mass parameters. 

In our simple minded analysis, we have used the partially coded field of view of the 
SPI detector given by Qpcfov ~ 35°. This corresponds to a solid angle Q = 27r(l — 
cos(8pcfov/2)) ~ 0-29 and we have used this value of the solid angle for our numerical 
calculation. The dark matter decay line flux has been calculated using various profiles 
mentioned in the previous section. In addition, we have used the minimal model described 
in Ref. 32[, for which the dark matter density is constant within r whereas it is given by 



the model A2 of Klypin et al in Ref. 



The results obtained in Ref. 



22j for r > r Q . 



32| using the SPI data set is applicable to any model of 



decaying dark matter. We have used the results shown in figure 11 of Ref.[32[], for the 3a 
restrictions on the lifetime of the dark matter particle as a function of the X-ray photon 
energy . Using the minimal model for the dark matter profile as mentioned above, we 
translate this restriction on the lifetime into a restriction on the decay line flux as a function 
of the axino mass. The dark matter decay line flux is also a function of the angular distance 
between the direction along the line of sight and the direction towards the Galactic center. 
We have used <p = 13° (the inner Galaxy) for the calculation of the flux restriction. 

In figure [5] we show the flux restriction using the red zigzag curve. In addition, in the 
same figure we show the line flux as a function of the axino mass coming from an axino dark 



matter, for different values of f a . These curves are drawn for the minimal profile of 32|. 



From the figure we see that for the astrophysically unconstrained range of the Peccei- 



1 We sincerely thank Oleg Ruchayskiy for providing us with the data file corresponding to figure 11 of their 
paper. 
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FIG. 5: Flux of photons from the Milky Way inner Galaxy (eft = 13°) for the minimal profile 
described in the text. The red zigzag line is the lower limit of the flux coming from the negative 
search of a dark matter decay line using the data from SPI spectrometer on board the INTEGRAL 
satellite. The shaded region is the window on the axino mass where the axino can act as a warm 
dark matter candidate. 



Quinn symmetry breaking scale, namely, f a > 10 9 GeV (25j, axino mass m„ > 50 keV could 
be excluded by SPI data. 

On the other hand, the results depend on the profile one chooses for the dark matter in 
the Milky Way. In particular, this dependence is more pronounced in the direction of the 
inner Galaxy. In figure [6] we plot the flux for a fixed value of f a for the isothermal profile 
and different models of the NFW profile (the values of the different parameters appearing 
in these profiles are taken from 32J) as well as the minimal profile of 32[. The shaded 
region and the zigzag line have the same meaning as before. From the figure it is obvious 
that there could be some noticeable variation in the flux depending on the profile one uses. 
Accordingly this gives us different upper bounds on the axino mass parameter. However, 
this variation becomes more and more negligible for some of the profiles if one looks at the 
direction of the Galactic anti-center. 



11 




10 



-2 





CM 




o 10 



-4 




Favored A-, or B-| 
Max. disk A 2 
Max. disk B 2 
Battaglia 



40 



60 



80 



100 120 140 160 180 200 



mg (keV) 



FIG. 6: Flux of photons from the Milky Way inner Galaxy (</> = 13°) for all profiles for f a = 10 9 
GeV. As before, the red zigzag line is the lower limit of the flux coming from the SPI spectrometer 
observations. The shaded region is the window on the axino mass where the axino can act as an 
warm dark matter candidate. 

V. SUMMARY, CONCLUSIONS AND OUTLOOK 

We have attempted to constrain the parameter space of a warm dark matter candidate 
from the non-observation of galactic X-ray line flux at the X-ray telescopes such as Chandra 
and SPI. We are considering an R-parity violating SUSY scenario, where the gravitino and 
the axino are the two likely candidates of warm dark matter. Due to the violation of 
lepton number, they can have some two body decays into a (monochromatic) photon and 
a neutrino, although their mean lifetimes are larger than the age of the universe. These 
photons correspond to X-ray lines that are potentially observable at X-ray telescopes. 

We have found that gravitino decays in the above channel yield too small a flux to lead 
to any useful constraint, due to the suppression of gravitino interaction by the Planck mass. 
The axino, on the other hand, couples with a strength inversely proportional to the Peccei- 
Quinn symmetry breaking scale (/ a ) , which can lie in the range 10 9-12 GeV, and is therefore 
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more promising from the viewpoint of X-ray satellite observations. Our analysis is based on 
the assumption that R-parity violation, triggered by bilinear terms in the superpotential, is 
of such a degree as to accommodate the pattern of neutrino masses and mixing suggested 
by the solar and atmospheric neutrino data. On this basis, we find that the limits from 
Chandra and SPI disallow a sizable part of the yet allowed axino parameter space defined 
by the axino mass (m„) and f a . In particular, axino masses in the range ~ 50 - 100 keV 
tend to be disallowed by SPI. The best results follow from a study of X-rays from the inner 
part of the Milky Way. For galaxy clusters such as the Virgo cluster a similar analysis has 
been performed for Chandra with similar conclusions as in the case of Milky Way. 

If one uses the exact expressions for neutrino masses, there could be some SUSY param- 
eter space dependence of the above results. In particular, there could be some cancellations 
in the expression for sneutrino VEV in some regions of the parameter space. However, the 
general conclusions, based on the limits on the extent of R-parity violation, are unlikely to 
differ drastically from what have been derived above. 

It should be noted that the limits have some dependence on the profile of the dark matter. 
The flux from the inner part of the galaxy can vary significantly depending on the profile. 
For some profiles like NFW favored models A\ or B%, NFW parameters of Battaglia et al, 
there could be constraints for f a = 10 9 GeV and m 5 > 42 keV. Even constraints can be 
obtained for f a = 10 10 GeV from SPI data. On the whole, for the lower range of allowed 
f a , namely, 10 9-10 GeV, axino masses above 42-60 keV (depending on f a ) can be ruled out 
from the SPI data on Milky Way. 

At the end, we would like to remark that the Chandra observatory is a high resolution 
experiment, but its energy range for X-rays is restricted to 8 keV. Similarly, for SPI the 
energy range is 20 keV to 7 MeV. There are other experiments; among them, XMM-Newton 
has less resolution but energy range up to 12 keV 33|. Beppo-SAXX can detect X-rays up 
to 15 keV 34j. Old experiments like HEAO-I, and HEAO-III can measure X-rays up to 
10 MeV. HEAO-II (Einstein) can measure X-rays up to 20 keV 35J. Japan's fifth X-ray 
satellite Suzaku can carry out high resolution spectroscopic observations in a wide energy 



range of 0.3 keV to 600 keV 



361 ] . A global analysis based on all of these experiments can 



lead to further insight into decaying warm dark matter candidates in an R-parity violating 
scenario. 
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